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ABSTRACT 

In the last stages of a black hole merger, the binary can experience a recoil due to asymmetric 
emission of gravitational radiation. Recent numerical relativity simulations suggest that the recoil 
velocity can be as high as a few thousands kilometers per second for particular configurations. We 
consider here the effect of a worst case scenario for orbital and phase configurations on the hierarchical 
evolution of the massive black hole (MBH) population. The orbital configuration and spin orientation 
in the plane is chosen to be the one yielding the highest possible kick. Masses and spin magnitudes 
are instead derived self-consistently from the MBH evolutionary models. If seeds form early, e.g. 
as remnants of the first stars, almost the totality of the first few generation of binaries are ejected. 
The fraction of lost binaries decreases at later times due to a combination of the binary mass ratio 
distribution becoming shallower, and the deepening of the hosts potential wells. If seeds form at later 
times, in more massive halos, the retention rate is much higher. We show that the gravitational recoil 
does not pose a threat to the evolution of the MBH population that we observe locally in cither case, 
although high mass seeds seem to be favored. The gravitational recoil is instead a real hazard for (i) 
MBHs in biased halos at high-redshift, where mergers are more common, and the potential wells still 
relatively shallow. Similarly, it is very challenging to retain (ii) MBHs merging in star clusters. 

Subject headings: cosmology: theory - black holes - galaxies: evolution - quasars: general 



1. INTRODUCTION 

Today, MBHs a re ubiquitous in the nuclei of nearby 
galaxies (see, e.g.. iRichstone et al]|1998l) . If MBHs were 
also common in the past, as implied by the large pop- 
ulation of quasars, and if their host galaxies experience 
multiple mergers during their lifetime, as expected in the 
currently favored cold dark matter hierarchical cosmolo- 
gies, then MBH mergers are expected to be relatively 
common during the cosmic history. 

The results of recent numerical relativity simulations 
question whether MBH mergers represent a threat to the 
evolution of the MBH population that we observe today. 
In the final phase of a black hole merger, the inspiral 
is driven by emission of gravitational radiation. Grav- 
itational waves carry, in general, a non-zero net linear 
momentum, which establishes a preferential direction for 
the propagation of the waves. As a consequence, the cen- 
ter of mass of the binary recoils in the opposite direction 
l|Redmount fc Ree"slll989f ). possibly causing the ejection 
of MBHs from the potential wells of their host galax- 
ies (e.g., IMadau et al.l 12004 iMadau fc QuataertJ 12004 
Merritt et all 12004 IVolonteri fc Reed 120061 : iHaimanl 
2004 ISchnittman fc Buonannoll2007l ). 

The calculations of the magnitude of the recoil 
have been oscillating in literature by a few orders of 
magnitude fro m the early ca lculations in the New- 
tonian regime jlFitchettl Il983f). to the latest analyti- 
cal approaches dFavata et al.ll2004t iBlanchet et al.ll2005t 
iDamour fc Gopakumarl |20*06T ). The advent of numeri- 
cal relativity is now leading to a convergence in the es- 
timates of the r ecoil. Schwarzsch ild, i.e., non-spinning, 
black holes (e.g. JBaker et al.ll2006l ) are expected to recoil 
with velocities below 200 kms -1 , and a similar range is 
expected for black holes with low spins, or with spins 

1 University of Michigan, Astronomy Department, Ann Arbor, 
MI, 48109 



(anti-) aligned with the orbital axis (Figure [T]). How- 
ever, when the spin vectors have opposite directions 
and are in the orbital plane, the reco il velocity can be 
as large as a few thousands kms" 1 (|Campanelli et al.l 
I2007at iGonzalez et alj[2007t ICampanelli et al.ll2007bf l. 

In this paper we assess if standard models for MBH 
evolution (mass and spin) can accommodate a high re- 
coil, and when such high recoils are indeed creating haz- 
ards during the hierarchical growth of the MBH popu- 
lation. We trace the mass and spin evolution of MBHs 
from early times, within a plausible scenario for the hi- 
erarchical assembly, growth, and dynamics of MBHs in a 
ACDM cosmology, that has been shown to capture many 
features of the MBH population (e.g., luminosity func- 
tion of quasars, MBH mass density, core formation due 
to MBH binary mergers). The main features of the mod - 
els have been dis c ussed i n IVolonteri et alJ (120 03. 2005); 
IVolonteri fc Reesl (|2006l ): ISesana et al l (|2007l ). and ref- 
erences therein. During a MBH merger, we assume the 
orbital and phase 2 configuration that leads to the largest 
recoil magnitude. Black hole masses and spins mag- 
nitudes are instead self-consistently determined by the 
merger and accretion history of the MBHs. 

We compare here two different models of MBH for- 
mation, one that predicts MBHs to form early, and 
with low mass, as remnants of the first generation of 
stars (PopIII stars, cfr. model VHM in Sesana et al. 
2007), one which predicts MBHs to form later and 
with larger masses (cfr. model BVRlf in Sesana et al. 
2007), in halos with lo w angular momentum, prone to 
large scale instabilities (jBegelman et al.l [20061 ) . We re- 

2 Campanclli et al. (2007a) also pointed out that the recoil ve- 
locity is not uniquely determined solely by the magnitudes and 
relative directions of the black hole spins in the orbital plane. The 
recoil was found to vary sinusoidally with the angle between the 
orbital plane (parallel to the black hole spins) and the initial linear 
momenta of the holes. 
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fer the reader to ISesana et al.l (|2007f ) for an exhaustive 
description of the models (including accretion proper- 
ties and dynamical evolution). We have assessed mod- 
els with MBH formation more/less efficient than the 
reference cases. In case of very efficient MBH forma- 
tion, the mass density on MBHs at z ~ 3 is much 
large r than expected from the observed quasar evolu- 
tion (|Yu fc Tremaind I2002t [Merlonil I2004D . In models 
w here MBHs form in very small numbers (e.g., BVRhf 
in lSesana et al.ll2007f ). mergers are very rare events, and 
happen at later times, with small mass ratios. Clearly, 
the gravitational rocket is much less of a hazard in 
models which start from fewer seeds at early times. 




q=( M MBH 2 / M MBHl) 

Fig. 1. — Recoil velocity of spinning black holes as a function 
of binary mass ratio, q. Solid curves: spin axis aligned (or anti- 
aligned) with the orbital angular momentum (Baker et al. 2007). 
Dashed curves: spin axis in the orbital plane (Campanelli et al. 
2007). For every mass ratio we plot the combination of spins, d\ 
and a,2 which minimizes (lower curves) or maximizes (upper curves) 
the recoil velocity. 

2. MERGING MBH MASS RATIOS AND SPINS 

During a galactic merger, dynamical friction drags in 
the satellite, along with its central MBH, towards the 
center of the more massive progenitor. The orbital evo- 
lution of the galaxy pair leads also to the destabilization 
of gas in the merging systems. In simulations, this gas is 
seen to accumulate towards the centers of the system and 
fuel the MBHs in both galaxies fe.g.. | Di Matteo et al.l 
[2005t iKazantzidis et al.ll200l iMaver et al.ll2006fl . When 
the orbital decay is efficient, the satellite hole moves to- 
wards the center of the more massive progenitor, leading 
to the formation of a bound MBH binary. 

The efficiency of dynamical friction decreases with the 
mass ratio of the merging galaxies, so nearly equal mass 
galaxy mergers ("major mergers") lead to efficient binary 
formation within short timescales (i.e., shorter than the 
Hubble time), while in "minor mergers" it can take longer 
than an Hubble time to drag the satellite hole to the cen- 
ter. This effects must be convolved with the mass ratio 
probability distribution. As the mass function of halos 
(and galaxies) is steep, the probability of halo mergers 
decreases with increasing mass ratio, that is, dynamically 
efficient major mergers are rare, while minor mergers, in- 
efficient in forming MBH binaries, are common. 

The combination of those effects leads to a distri- 



bution of mass ratios for merging black holes, which 
is at most redshifts dominated by systems with q = 
Mmbh2/Mmbhi < 1 in the decade around q = 0.1. 
At early times, equal mass MBH mergers are more com- 
mon, as most mergers involve MBHs which had no time 
to accrete much mass yet. At later times binaries with 
small mass ratios are instead much more common, as 
some MBHs that started their orbital decay during high 
rcdshift minor mergers (with an orbital decay timescale 
comparable to the Hubble time) finally find their way to 
the central MBH. The lower redshift mas s ratio distri- 
butio ns are therefore much shallower (cfr. ISesana et al.l 
I2007T ) : high-speed recoils are common at high redshift, 
but the typical strength of the re coil decreases with t ime- 
Note that the recoil estim ates by iBaker et al.l (|2007f) and 
I Campanelli et al.l (|2007al fbl) both imply velocities below 
a few kms" 1 for mass ratios q < 0.01. 

The other two variables determining the recoil veloc- 
ity (as we have fixed the relative inclination of the spins 
with respect to each other and the orbital plane) are 
the MBHs spin magnitudes 3 . MBH spins, in their cos- 
mic evolution, are determined by the sequence of MBH 
mergers and accretion episodes that build up the MBH 
mass. In general, mergers of similar mass MBHs in- 
crease the spin, while the capture of smaller compan- 
ions in randomly-oriented or bits tends to spin holes down 
(|Hughes fc Blandfordl 12003). Given our distribution of 
MBH binary mass ratios, binary coalescences alone do 
not lead to a dominant s pin-up or spin-down of MBHs 
(cfr. IVolonteri et al.l l2005h . Gas accretion, on the other 
hand, seems to be more efficient in determining the distri- 
bution of MBH spins, due to the alignment between spin 
and angular momentum of the accretion disc induced by 
the Lense-Thirring precession. 

In our simulations, we assume that in an accretion 
episode the MBH spin is initially misaligned with respect 
to the angular momentum of the disc by a random angle, 
and that the alignment i s relatively inefficient, following 
the scheme described in IVolonteri et al.l (|2005l ) (eq. 17- 
19; see Volonteri, Sikora & Lasota 2007 for a more com- 
prehensive model). The spin distributions that we derive 
are dominated by Kerr MBHs 4 , and are steeper for the 
spin of the primary (i.e. the most massive hole in the 
binary), as it typically grows by a larger factor during a 
given merger. 

3. DISCUSSION 

3.1. Signatures on the local MBHs population 

Our set of simulations allows to assess the influence of 
the gravitational recoil under realistic, although rather 
pessimistic, assumptions. The fraction of merging MBHs 
that are actually ejected from their hosts, i.e., that have 
a recoil velocity larger than the escape speed from the 
host, has a strong dependence on redshift and on the 
MBH seeds formation mechanism (Figure[2]). The escape 
velocity calculation assumes an NFW halo for the DM 

3 We express the spin through the dimensionless parameter a, 
S = aGm\ K /c, < a < 1 

4 For simplicity, we will assume in the following that a MBH may 
be spun up to a maximum equilibrium value of a = 0.998 (Thorne 
1974), although we note that relativistic magnetohydrodynamic 
simulations for a series of thick accretion disc models converge at 
a spin equilibrium of d s; 0.93 (Gammie, Shapiro, & McKinney 
2004). 
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component. The baryonic component is modeled as an 
isothermal sphere truncated at the radius of the MBH 
sphere of influence, that is ri n f — GMmbh/^q, where 
Mmbh is the MBH mass and a g is the velocity dispersion. 

If seeds form early, in small proto-galaxies, almost the 
totality of the first few generation of binaries are ejected 
into the intergalactic medium. The fraction of "lost" 
binaries decreases at later times due to a combination of 

(i) the mass ratio distribution becoming shallower, and, 

(ii) the hierarchical growth of the hosts. The ejected 
fraction decreases below 50% at z ~ 5. If seeds form at 
later times, in more massive halos, the retention rate is 
much higher from the beginning, and only about 40-50% 
of binaries are ejected at all redshifts higher than ~3. 
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Fig. 2. — Merger rate of binaries (solid histogram) and ejection 
rate of these binaries (dashed histogram), a binary is defined ejected 
if the recoil velocity is larger than the escape velocity frm the host. 
Upper panel: popIII remnants; lower panel: direct collapse. 

Let us now quantify if these ejection rates indeed dis- 
turb the global evolution of the MBH population. The 
number of MBHs in binaries is a very small fraction of 
the population at most times (Figure [3]) in all the con- 
sidered cases. In fact, only about 40% of MBHs ever 
experience a merger, and MBHs evolve in isolation for a 
large part of their lifetime. As, especially at high red- 
shift, binaries represent the exception, rather than the 
rule, the possible ejection of most binaries before z ~ 5 
is not a threat to the evolution of the MBH population 
that has been detected in nearby galaxies. Provided that 
one o r few MBHs per galaxy merger tree (M enou et alJ 
l200ll ) survive until z ~ 5, Soltan-type arguments, which 
compare the mass density of MBHs in today quiescent 
galaxies with the mass density inferred from luminous 
quasars, imply that accretion alone can grow the MBH 
up to the supermassive variety (Yu & Tremaine 2003). 

Figure 2] compares the simulated MBHs that are left 
in galaxies at z — to the current observational sample 
(samples from Ferrarese & Ford 2005; and Tremaine et 
al. 2002, augmented for Lauer et al. 2006 additions; see 
Valluri et al. 2004 for uncertainties in the MBH mass 
measurements). The simulated sample broadly follows 



the same scaling with galaxy velocity dispersion, a g , of 
the observational sample. By visual inspection, the sim- 
ulated MBHs display some "outliers" with masses below 
the predictions of the MBH mass - a g scaling. Given 
the small observational sample, and the uncertainties 
on the MBH masses, we ran a Montecarlo test to com- 
pare the obse rvational datasets to our simulations (cfr. 
IWvithd 120061 ): we conclude that statistically we can- 
not rule out that the simulated sample is compatible 
with observations (at the 2-a level), at least for MBHs 
hosted in galaxies with velocity dispersion larger than > 
100 km s . Additionally, observations are biased against 
the detection of "under-massive black holes" , with re- 
spect to the expectations of the MBH mass - a g scaling, 
given the scaling of the radius of the sphere of influence. 
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Fig. 3. — Number density of MBHs (upper curves) and binaries 
(lower curves). The fraction of MBHs in binaries is very small at 
most times. Upper panel: direct collapse. Lower panel: popIII 
remnants. 

Models with large MBH seed masses appear to be 
favored in a high recoil velocity scenario. However, 
note that the l-cr di s persions of the two samples, 
(jFerrarese fc Fordl 120051 : iTremaine et"aH l2002f ) are com- 
patible, but not identical. The fraction of outliers can 
differ up to a factor of two between the two samples for 
galaxies with velocity dispersion below 200 km s -1 . The 
uncertainties and small size of the observational sam- 
ple contribute to prevent to draw firm conclusions: a 
larger sample of measured MBH masses, which can be 
available to future instruments, such as Giant Magellan 
Telescope, can help assess the relative importance of the 
gravitational rocket. 

We recall here that we have adopted the orbital con- 
figuration yielding the highest recoi l veloc ities, which, 
as pointed out by iBogdanovic et ail (120071) is prob ably 
rather uncommon. Also. ICampanelli et al.l (|2007aD dis- 
cuss how the sinusoidal dependence on the orbital phase 
implies a reduction of y2 in the root-mean-square recoil 
velocity, in case of optimal orbital configuration with ran- 
dom phase configuration. The fraction of outliers that we 
find represents therefore an upper limit, and our results 
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emphasize also the necessity of further investigations on 
the orbital and phase configuration expected in merging 
MBH binaries. 
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Fig. 4. — The Msmbh - velocity dispersion (a c ) relation at z = 0. 
Every circle represents the central MBH in a halo of given tr c . 
Observational data are marked by their quoted errorbars, both in 
<t c , and in Msmbh (Tremaine et al. 2002). Left: popIII remnants; 
right: direct collapse. Upper panels: fraction of galaxies at a given 
velocity dispersion which do not have a central SMBH. 

3.2. Special cases 

While Figure [2] is reassuring with respect to the "av- 
erage" MBH, as the simulations sample a statistical en- 
semble of galaxy halos, it is not representative of very 
rare, highly biased regions, where MB H forma t ion is 
more common, and so are binary MBHs. iHaimanl (pOOi ) 
pointed out that the recoil can be indeed threatening 
the growth of the MBHs that are believed to be power- 
ing the lumi nous quasars at z — 6 detected in the Sloan 
survey (e.g.. iFan et al.1 120011 ). In fact, in such a biased 
volume, the density of halos where MBH formation can 
be efficient (either by direct collapse, or via PopIII stars) 
is highly enhanced. The net result is an higher fraction 
of binary systems, and binarity is especially common for 



the central galaxy of the main halo. While the "average" 
MBH experiences at most one merger in its lifetime, a 
MBH hosted in a rare exceptionally massive halo can ex- 
perience up to a few tens mergers, and the probability of 
ejecting th e central MBH, ha l ting i ts growth, is 50-80% 
at z > 6. IVolonteri fe Reesl (|2006h demonstrated that 
this result is not dependent on the "optimal" orbital and 
phase configuration. This implies that MBHs at high 
redshift did not mainly grow via mergers, as proposed 
in works that did not take the gravitational recoil int o 
consideration (e.g.. lPelupessv et al.ll2007tlLi et al. ll2006). 

Another special case, at the other very end of the MBH 
mass scale, is represented by merg ing MBHs in star clus- 
ters. It has been propose d (e.g., iPortegies Zwart et al.1 
I2004t iGiirkan et atl |2004[ ) that the merging of main- 
sequence stars via direct physical collisions can enter 
into a runaway phase, and form a very massive star, 
which may then collapse to form a MBH. A more ex- 
otic outcome of the process can happen in star clusters 
with a relatively large (but well within the expected val- 
ues) pr imordial binary fra ction: the formation of two 
MBHs (jGurkan et al.H2006l ). Star collisions likely trans- 
fer angular momentum to the very massive star, and 
a rotating very massive star pla usibly collapses into a 
spinning MBH (|Saiio et al J 120021 ). If the two star clus- 
ter MBHs form with a mass ratio larger than q = 0.1 , 
when they proceed to coa lescence (jGultekin et al.ll2~004 
20061; iFregeau et all [2006). the magnitude of the recoil 
for most configurations is larger than 50 kms -1 , lead- 
ing to a high escape probability. If formation of binary 
MBHs is common in star clusters, their ejection is highly 
likely as well. A similar fate would await MBHs bina- 
ries forming follow ing the merger of two star cl uster, 
each with a MBH (|Amaro-Seoane fc Freitagll2006[ K Al- 
though in both cases the MBH binary would be a source 
of gravitational radiation, which might allow to confirm 
the predictions of the previous scenario for MBH for- 
mation in cluster, after the ejection the MBHs would 
join the ve ry difficult to detect population of wa nder- 
ing MBHs (Volont erTfc Pernall2005l: IM apelli et all 120061 : 
iKuranov et afll2007f ). (|Volonteri et al.ll2007| ) 
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